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Theorie mathematique des modMes de MHD reduite pour les 

plasmas de fusion 

Resume : L'etablissement de modeles de MHD reduite est formule comme un exemple de la 
theorie generale des limites singulieres des systemes h 5 rperboliques. Cette formulation permet 
d'utiliser les resultats generaux de cette theorie et de prouver rigoureusement que les modeles 
de MHD reduite sont une approximation valide du modele complet. En particulier, la conver¬ 
gence des solutions du modele complet vers les solutions d'un systeme reduit est demontree. 

Mots-cles : Analyse asymptotique, systemes h 5 rperboliques, limite singuliere, MHD, Plasmas 
de fusion, Tokamaks 
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1 Introduction 

Magnetohydrod 5 mamics (MHD) is a macroscopic theory describing electrically conducting flu¬ 
ids. It addresses laboratory as well as astrophysical plasmas and therefore is exfensively used 
in very differenf confexfs. One of fhese confexfs concerns fhe sfudy of fusion plasmas in foka- 
mak machines. A fokamak is a toroidal device in which hydrogen isotopes in fhe form of a 
plasma reaching a femperafure of fhe order of fhe hundred of millions of Kelvins is confined 
fhanks to a very sfrong applied magnetic field. Tokamaks are used fo sfudy confrolled fusion 
and are considered as one of fhe mosf promising concepfs fo produce fusion energy in fhe near 
fufure. However a hof plasma as fhe one presenf in a fokamak is subjecf fo a very large num¬ 
ber of insfabilifies fhaf can lead fo fhe end of fhe existence of fhe plasma. An imporfanf goal 
of MHD sfudies in fokamaks is fherefore fo defermine fhe sfabilify domain fhaf consfrainfs fhe 
operational range of fhe machines. A secondary goal of fhese sfudies is fo evaluate fhe conse¬ 
quences of fhese possible insfabilifies in term of heaf loads and sfresses on fhe plasma facing 
componenfs of fhe machines. Numerical simulations using fhe MHD models are fherefore of 
uffermosf imporfance in fhis field and fherefore fhe design of MHD models and of models be¬ 
yond fhe sfandard one (e.g incorporafing fwo-fluid or kinefic effecfs) is fhe subjecf of an infense 
acfivify. 

However, fhe MHD model is a very complex one : if confains 8 independenf variables, fhree ve- 
locify componenfs, fhree componenfs of fhe magnefic field, densify and pressure. Alfhough fhe 
sysfem is h 5 rperbolic, if is nof sfricfly hyperbolic leading fo fhe existence of possible resonance 
befween waves of differenf f 5 rpes and moreover fhe MHD sysfem has fhe addifional complexify 
of being endowed wifh an involufion. An involufion in fhe sense of conservafion law sysfems is 
an addifional equation fhaf if satisfied af f = 0 is safisfied for all t >0 13. For all fhese reasons, 
approximafions and simplified models have been designed bofh for fheorefical sfudies as well 
as numerical ones. In fhe field of fusion plasmas, fhese models are denoted as reduced MHD 
model^ These models initially proposed in fhe 70' ETIH have been progressively refined fo in¬ 
clude more and more physical effecfs and correcfions Il28l l3l l22l[12l . In particular, some earlier 
models conserve a non-sfandard energy and in some modeling works, special affenfion have 
been paid fo insure fhe conservafion of fhe usual energy e.g IIT^ ItI ITtII (see also 111). Af presenf 
fhe liferafure on fhe physics of fusion plasma concerned by reduced MHD models is huge and 
confains several hundred of references. From a numerical poinf of view, several well-known 
numerical codes (e.g Il211l , JH ) used routinely for fusion plasma sfudies are based on fhese re¬ 
duced models. Acfually, while fhere is a definife fendency in fhe fusion plasma communify fo 
use full MHD models e.g 13, 1T3, ITOl , a large majorify of non-linear simulations of fokamak 
plasmas have been and still are conducted wifh fhese approximafions. 

Until recenfly, reduced MHD models have nof affracfed a lof of inferesf in fhe mafhemafical 
or numerical analysis liferafure. One can cife 13 and || 8 | fhaf have shown fhaf fhese models 
can be interpreted as some special case of "Galerkin" mefhods where fhe velocify and magnefic 
fields are consfrained fo belong fo some lower dimensional space. This inferprefafion is also 
implicif in fhe design of fhe M3D-C1 code 113] where insfead of fhe usual projection on fhe 
coordinate sysfem axis, fhe equafions governing fhe scalar componenfs of fhe vecfor fields are 
obfained by special projections fhaf allow fo recover reduced models. 

In fhis work, we adopf fhe differenf poinf of view of asympfofic analysis and show fhaf 
reduced MHD models can be undersfood as a special insfance of fhe general fheory of singular 
limif of hyperbolic sysfem of PDFs wifh large consfanf operafors. This formulation allows fo use 

^while the standard MHD model is by contrast designated as the full MHD model 
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the general results of this theory and to prove rigorously the validity of fhese approximations 
of fhe MHD equations. In particular, if is proven here, we believe for fhe firsf time, fhaf fhe 
solufions of fhe full MHD sysfem converge fo fhe solufions of an appropriafe reduced model. 

This paper is organized as follows : Firsf, we recall fhe general fheory of singular limifs of 
quasi-linear h 5 rperbolic sysfem wifh a large paramefer. In fhe third section, we show how this 
general framework can be used fo analyze reduced MHD models. Finally, we conclude by some 
remarks on possible exfensions of fhe presenf work. 


2 Singular limit of hyperbolic PDEs 

2.1 General framework 

In this section, we are concerned with the behavior when e —^ 0 of fhe solufions of h 5 rperbolic 
sysfem of PDEs of fhe following form : 

I Ao{W,e)dtW + j:^A,{W,e)d,^W = 0 

\ W{0,x,e) = Wo{x,e) 

Here VF G 5 c is a vector function depending oi = 1, • • • , d) where d is fhe space 

dimension while fhe ^O; Aj, Cj are square N x N mafrices. Due fo fhe presence of fhe large 
coefficienf 1 /e multiplying fhe operafor {■), we may expecf fhe velocify of some waves 

present in Q to become infinite and therefore, for a solufion fo exisf on a 0(1) fime scale, if 
has fo be close in some sense fo fhe kernel K = {W G IR^s.tJ^j CjdxjW = 0} of fhe large 
operafor. The limit system obtained from Q is fherefore a singular limit since the constraint 
W G K may change the h 5 rperbolic nature of fhe sysfem Q. A profofypical example of fhis be¬ 
havior is given by fhe incompressible limif of fhe hyperbolic equations governing compressible 
Euler flows where fhe propagafion af infinite speed of fhe acousfic waves gives rise fo an elliptic 
equation on fhe pressure coming from fhe global consfrainf V • tt = 0. 

The nafure of fhe singular limif depends on fhe initial dafa. Using fhe ferminology of Scho- 
chef ||25]| fhe limif is called "slow" if the initial data makes the first time derivatives at time t — 0 
stay bounded as e —0. The term "well-prepared initial data" is also used to qualify fhis sifua- 
tion. In fhis case, under appropriafe assumpfions, fhe solufions exisf for a fime T independenf 
of e and converge fo fhe solufions of a limif sysfem when e —> 0. 

In fhe opposife case, denofed as a "fasf" singular limit, dtW is not 0{1) at time zero and 
fast oscillations developing on a 1/e time scale can persist on the long time scale. Solutions of 
fasf singular limif cannof converge as e —>^ 0 in fhe usual sense since fhe fime derivative of fhe 
solufion is of order 1/e. In fhis case, convergence means fhe exisfence of an "averaged" limif 
profile W(f, r, x) such fhaf W(t, x, e) — W{t, t/e, x) ^ 0 wifh e. The question of fhe exisfence 
of fasf singular limif is in parficular examined in 125] . A review article summarizing resulfs on 
fhis subjecf wifh a special emphasis on fhe low Mach number limif is m. 

In fhis work, we will be mainly concerned by fhe slow case. Even in fhis case, fhe exisfence 
for a fime independenf of e and fhe convergence of fhe solufions fo fhe solufions of a limif 
sysfem may require additional assumptions on the structure of Q. Beginning wifh fhe earlier 
works in the 80' of Klainerman and Majda ITSl 1141 [T9l and fhose of Kreiss and his co-workers 
1T^ |21, fhese quesfions have been examined in several works 1^124] with the main objective to 
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justify the passage to the incompressible limit in low Mach number compressible flows. Several 
extensions of these works for viscous flows or general hyperbolic-parabolic systems are also 
available. Again we can refer to ||T| for a review. 

The following theorem (see IIT9]l , chapter 2) summarizes the main results of these works in a 
form suitable for our purposes : 


Theorem 11. Assume that: 

1. Conditions on the initial data : Wo(x, e) = Wq{x) + eWQ{x, e) 

(a) Wq{x) and Wl{x,e) are in 

(b) E, C,d,W° = 0 

(c) \\Wl{x,e)\\s<Cte 
1. Structure of the system 

(a) The matrices Aq, Aj and Cj are symmetric 

(b) Aq is positive definite at least in a neighborhood of the initial data 

(c) Aq and Aj are C® continuous for some s > [n/2] -|- 2, where n is the number of spatial 
dimensions 

(d) The Cj are constant matrices 

(e) The matrix Aq{W, e) = Ao{eW) 

then the solution W{t,x, e) of system Q with the initial data satisfying condition 1 is unique and exists 
for a time T independent of e. In addition the solutions W{t, x, e) satisfy : 

\\W{t,x,e) -W°{t,x)\\,_i <CeiortG [0,T] 

where W°{t, x) is the solution of the reduced system : 

( Ao(0)9tTr° + E, = 0 

{ ( 2 ) 

[ W°{0,x) = W°{x) 

Proof. : The proof of this result can be found in 119] , chapter 2, Theorems 2.3 and 2.4. We do not 
repeat this proof here but briefly comment on some of their aspects: The assumptions 2.(a) et 
2.(b) simply means that system Q is a quasi-linear symmetric hyperbolic system in the sense of 
Friedrichs. The uniqueness and existence of solution on a finite time T > 0 can then be estab¬ 
lished by classical iteration techniques relying on energy estimates (see for instance HSl or lT9l l. 
However the presence of the large coefficient 1 je could possibly make this time of existence e- 
dependent and shrinking to 0 with e. Assumption 2.(d) ensures that this will not be the case 
since the matrices Cj being constant, the large terms will not contribute to the energy estimates. 

The assumption 2.(e) Aq = A^^eW) allows to bound its time derivative independently of e : 
Since we have 

dMo{eW) = ^ed,W = -^ eAf^[A,d,W + ^d^C^W] 

The £ and 1 /e terms balance together and give an estimate independent of e. 
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The assumptions l.(b) and l.(c) means that the initial condition is sufficiently close to the 
kernel of fhe large operator to ensure fhaf fhe time derivative dtW (0, x, e) is bounded in 
independenfly of e. This condition implies fhaf fhe inifial dafa are "well-prepared" and will nof 
generafe fasf oscillations on a 1/e time scale. □ 

2.2 Reduced limit system 

Even if ||^ provides a complete description of fhe behavior of fhe solutions of fhe original sys¬ 
tem as £ fends to 0, fhe limif system confains as many unknowns as fhe original one. Acfually, 
one may even consider fhaf if confains more unknowns as fhe firsf order correcfion have 
also fo be computed. In practice, fhis largely depends on fhe specific sysfem considered as some 
lines of fhe mafrices Aj 0) may be identically zero and/or fhe evaluation of some ferms of 
fhe firsf order correcfion can be complefely obvious. However, if can be inferesfing fo derive 
from ||^ a "reduced" sef of equations confaining less unknowns by eliminating fhe firsf order 
correcfion W^. A particularly pleasanf framework fo consfrucf such a reduced sysfem is fhe 
following : 

Assume fhaf fhe kernel K = {W € IR^CjdxjW = 0} have dimension n < N and 
can be paramefrized by a linear operator wifh consfanf coefficienfs such fhaf: 

VWgKcM^, W = M{oj) 

Since fhe operafor L = have consfanf coefficienfs, A4{uj) is also a differential op- 

erafor of order 1 wifh consfanf coefficienfs fhaf can be written : 

d 

= +PoV (3) 

where fhe mafrices {Pj]j = 0,c?} are reef angular N x n consfanf mafrices. Then consider fhe 
adjoinf operafor Ai* from fo IR” safisfying 

(Af(w), W) = {uj,M*W) 

The operafor A4* is an "armhilafor" for fhe L operafor in fhe sense fhaf 

M*'L = 0 


Indeed we have: 

= (LW,Mu;) = -(Vte,LMa;) = 0 

since fhe Cj being symmefric mafrices, L is a skew-symmefric operafor. 

From ||^ Af * has fhe explicif expression : 

d 

M*iW) = -J2 W + PlW ( 4 ) 

where Pj; j = 0, • • • ,d are reef angular n x N mafrices, fransposes of fhe Pj. 

Wifh fhe operafors Ai and A = Ai* ai hand, a reduced sysfem of equations can be obfained 
by left multiplying 0 by fhe armhilafor A for functions W = AHui). In fhis operafion, fhe 
firsf-order correcfion Cj^xj vanishes and we obfain wifh W = Af (uj) : 

f AAQ{0)Mdt(^ + J2j-^^ji-^{‘^):0)Mdx^uj = 0 

\ a;(0, x) = uj°{x) 
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that is an autonomous system for the reduced variable uj € fR”. Note that to obtain ||^, we 
have used the fact that At being a linear differential operator defined by consfanf mafrices Pj, 
if commufes wifh fhe time and spafial derivatives. 

Nofe also fhaf spafial derivatives are "hidden" in fhe definifion of fhe operafors A and A4. 
Therefore in confrasf wifh fhe equafions (|^ fhaf is a firsf-order differential sysfem, ||^ defines 
a third-order differential sysfem of equafions (see sectio n |3. 2.3 for fhe concrefe example of re¬ 
duced MHD sysfem). The choice of using (|^ insfead of (|^ as a basis for a numerical mefhod is 
fherefore problem dependenf and in pracfice l|^ can be more difficulf fo approximafe fhan fhe 
original limif sysfem. 


3 Application to reduced MHD 

3.1 The ideal MHD system 

We now proceed fo show how fhis general framework can be applied fo fhe MHD equafions and 
begin fo recall some basic facfs abouf fhis sysfem. In fhe sequel, we will make fhe assumption 
fhaf fhe flow is barofropic, fhaf is fhe pressure is only a function of fhe densify. This assumpfions 
includes isenfropic as well as isofhermal flows. 

The ideal MHD sysfem can be wriffen under many differenf forms. Since fhe general fheory we 
have described make use of fhe symmefry of fhe jacobian mafrices, we use here a symmefric 
form of fhe sysfem : 


p^u + V{p + B^/2)-{B.V)B 

= 0 

( 6 . 1 ) 

— B — {B.\7)u B\7.u 

= 0 

( 6 . 2 ) 

1 D 

— 77713 + V.M 

7P JJt 

= 0 

(6.3) 


In fhese equation, u is fhe velocify, B fhe magnetic field and p is fhe pressure. The densify p is 
relafed fo fhe pressure by a sfafe law p = p{p), for insfance fhe perfecf gas sfafe law fhaf wrifes 
p = A{p /where A and 7 are consfanf and s is fhe (here consfanf) enfropy. The nofafion 
D./Dt sfands for fhe maferial derivative fhaf is defined by D ■ jDt = dt ■ -l-(tt.V)-. 

To sysfem ||^ one musf add fhe involution : 

V • B = 0 (7) 

and if is easily checked fhaf if 10 is verified af f = 0 , if is verified for all f > 0 . 

The sysfem 0 is h 5 rperbolic, ifs Jacobian has real eigenvalues and a complefe sef of eigen¬ 
vectors. However, if is nof a sfricfly h 5 rperbolic sysfem since some eigenvalues may coincide. 
Aparf from waves moving wifh fhe maferial velocify, if is usual fo splif fhe sef of MHD eigen¬ 
values and associated waves info fhree groups, fhaf are defined as : 
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Fast Magnetosonic waves : 


A^ = u.n ± Cf 

with 

II 

liVt^ + vl + ) 

(8.1) 

Alfen waves : 





A^ = u.n ± Ca 

with 

II 

(B.nf/p 

(8.2) 

Slow Magnetosonic waves : 
A^ = u.n ± Cs 

with 

_ 

Og — 

+ v\- 

(8.3) 


where va and Vt are defined hy : v\ = l-Bp/p and = yp/p. 


The velocity of these waves is ordered as follows : 

A| < Ai < A?. 


Fast and slow Magnetosonic waves are the equivalent of acoustic waves in fluid d 5 mamics. 
Alfen waves (sometimes also called shear Alfen waves) are of a different nature: The expression 
l|^2) shows that they do not propagate in the direction orthogonal to the the magnetic field. 
Actually in the direction orthogonal to the magnetic field, the speed of propagation of Alfen 
and slow magnetosonic waves is zero (in a frame moving with the material velocity) and only 
the fast magnetosonic waves survive. 


3.2 Large aspect ratio theory 

3.2.1 Geometry and coordinate system 

In this section, we are concerned with the model of the "straight tokamak" that consists of 
a slender torus characterized by a small aspect ratio e = q/Rq (see figure]^. In this model, 
the torus is approximated by a periodic cylinder of length 27ri?o arid of section of radius a. 
Some of the d 5 mamical effects that occur in a tokamak are well represented in this way and this 
model have been extensively used in theoretical studies to understand tokamak d 5 mamics. In 
particular, it is the model considered in 1271 to derive his original reduced model. 



Figure 1: Straight tokamak model: the slender torus is unfold to form a periodic cylinder 
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Now, let (^, 77 , Q be the usual cartesian coordinate system and let us introduce the normal¬ 
ized variables : 

( x = C/a 

S y = vla 

[ z = C/Ro 

In a way consistent with the underlying physical problem, the z direction will be denoted as the 
toroidal direction while the planes {x, y) are the poloidal sections. Note also that the introduc¬ 
tion of the normalized coordinates (x, y, z) corresponds actually to a two scale analysis : z the 
toroidal coordinate is scaled with Rq while the poloidal coordinates (x, y) are scaled with the 
small radius a. 

With these normalized coordinates, the expression of fhe spafial operafors becomes : 


aV/ 

df ^ d/ 

“ + 


dvz 

aV • V 

= Vj. • iix -f 

oz 

aV X V 

= (ez • V_L X Uj_)ez 


I' 


dvy dvx 
dz dz 


(9.1) 

(9.2) 


= - dxVz^y + {dxVy - dyVx)e2 + £(-d^Uj^ex -I- d^u^^ey) (9.3) 

wifh fhe defmifions : 


V = v±+ 


v± = 

Vx • uj_ = 


dvx 

dx 


dvy 

dy 


3.2.2 Scaling 

We now proceed fo scale fhe unknown variables. To recasf fhe equafions info an useful form, fhe 
usual procedure is fo wrife fhem in dimensionless form by scaling every variable by a characfer- 
isfic value. Here, in an equivalenf manner, we will consider fhe following change of variables: 


Magnefic field: 

B — — Gz + Bp — -So(ez H- sB) 

H 

( 10 . 1 ) 

Pressure : 

p = Po{p-\-£q) 

( 10 . 2 ) 

Velocify: 

U = £VaV 

(10.3) 

Time: 

I <.0 

II 

(10.4) 


where Bq is fhe reference value of fhe foroidal magnefic field on fhe magnefic axis {R = Rq) and 
p is a consfanf. In fhese expressions, va is the Alfen speed defined by v\ = po where po is 
some reference densify (for insfance, a characferisfic value of fhe densify on the magnetic axis). 
We choose for simplicity Pq = pqv\ (this only affects the value of fhe consfanf p). The imporfanf 
assumptions made in (TO) are : 

i) that the toroidal magnetic field dominafes fhe flow and fhat fhe poloidal field is of order e 
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with respect to the toroidal field : B = Bt + eBp. In tokamaks, the toroidal field is mainly due 

F 

fo exfemal coils and if varies f 5 rpically as Bp = where F is approximafely a consfanf and 

R 

R is fhe disfance fo fhe rofafion axis of fhe forus. In fhe model of fhe "sfraighf fokamak" and in 
fhe limif of small aspecf rafio q/Rq, fhis leads fo fhe following expansion of fhe magnefic field : 


B = -e^ 
R 


eBp = 


F-Fn 


Ro{l + ex) Ro{l + ex) 


Bz + eBp — BqI^Bz + eB) 


where Bq = Fq/Rq is fhe value of fhe toroidal magnefic field on fhe magnefic axis. Nofe fhaf 
B confains a foroidal componenf. This componenf is assumed fo be of fhe same order fhan fhe 
magnefic poloidal field. 


ii) fhaf fhe pressure flucfuafions are also of order e. Since fhe poloidal magnefic field is of 
order e wifh respecf fo fhe foroidal field, fhis means fhaf fhe poloidal plasma f3 paramefer is of 
order 1. In fhe physical liferafure, fhis sifuafion is referred fo as a "high" /? ordering 128] . 


iii) fhaf fhe velocifies are small wifh respecf fo fhe Alfen speed. Sfricfly speaking fhis assump- 
fion needs only fo be done for fhe perpendicular velocify. We adopf if for fhe full velocify vector 
in order fo simplify fhe presenfafion. 


we 


iv) fhaf we are inferesfed in fhe long fime behavior. Acfually, fhe assumpfion (TO) means fhaf 
are inferesfed in fhe long fime behavior of fhe sysfem wifh respecf fo fhe Alfen fime a/vA thaf 
represenf fhe fypical fime for a magnefosonic wave fo cross fhe fokamak secfion. (see section 


3.2.4 for some remarks on fhe shorf fime behavior of fhe sysfem on fhe fasf scale u/va)- 


Infroducing fhese expression info fhe MHD sysfem, we gef: 


P{P + £Q)] [^^ + + d,(g + B,)e^ + WpBy2 - d,B - {Bp ■ Vp)B 


+e{pv^dzV + dz{B'^ I2)b^ - B^d^B) 


-Vp{qFB,) 

e 


= 0 


^B± + {vp ■ \'p)Bp - {Bp ■ \'p)vp + BpVp ■ vp - d^vp 

OT 

+e{v^d^Bp - B^dzVp + d^v^^Bp) = 0 


( 11 . 2 ) 


( 11 . 3 ) 


A, 

dr 


1 


i{P + sq) 

dr 



= 0 


( 11 . 4 ) 

( 11 . 5 ) 
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If one introduces the variable W = {vj;,Vy, Vz,Bx, By,Bz,qY, the previous system can be 
written as 

A^{eW)drW + Y, eW)d^^W + - Y = 0 ( 12 ) 


where the matrices Aq, Aj {W, eW) are defined by : 


— 


/ pVx 
0 
0 
0 
By 
Bz 



( Ph 

O3 

0 

\ 


Aq = 

O3 

h 

0 

1 




1 

O3 

'r{p-\-eq) 7 


0 

0 

0 

By 

Bz 

0 

pVx 

0 

0 

-Bx 

0 

0 

0 

PVx 

0 

0 

-Bx 

0 

0 

0 

Vx 

0 

0 

0 

-Bx 

0 

0 

Vx 

0 

0 

0 

-Bx 

0 

0 

Vx 

0 


A, = 


V 


epvz 

0 

0 

-(l + eBz) 
0 
0 
0 


V 

0 

0 

0 

/ 

pVy 

0 

0 


0 

pvy 

0 


0 

0 

pvy 


-By 

Bx 

0 


0 

0 

0 


0 

Bz 

-By 

V 

0 

0 

0 


0 


0 


epvz 


0 


0 


epVz 


0 


eBx 


-By 

Bx 

0 

Vy 

0 

0 

0 


0 

0 

0 

0 

Vy 

0 

0 


0 

0 

Bz 

-By 

0 

0 

Vy 


-(1 + cBz) 
0 


SBy 

0 

1 


-(l + eBz) 
0 

sBx 

EVz 

0 

0 

0 


0 

-(l + sBz) 
SBy 
0 

evz 

0 

0 


')(p + £q) / 

0 \ 
0 
0 
0 
0 
0 

'•’y 

'yip + sq) ] 
0 


0 

0 

0 

0 

evz 

0 


r- 

'^{p + eq) / 


while the constant matrices Cj are given by : 


0 0 0 0 0 1 1 


Cx = 


0 0 0 0 0 0 


/ 

0 

0 

0 

0 

0 

0 

0 



0 

0 

0 

0 

0 

1 

1 


0 

0 

0 

0 

0 

0 

0 



0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 


Cy — 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 



0 

0 

0 

0 

0 

0 

0 


1 

0 

0 

0 

0 

0 

0 



0 

1 

0 

0 

0 

0 

0 

V 

1 

0 

0 

0 

0 

0 

0 

/ 

V 0 

1 

0 

0 

0 

0 

0 / 


This form makes apparent that the ideal MHD system can be put under the general form stud¬ 
ied in section 2.1 Therefore, the general results obtained in this section can be applied and we 
have 


Theorem 2 1. Assume that the initial velocity, magnetic field and pressure are defined by : 

( u{0,x)/Va = + ev^{e,x)) 

< B{f,x)/BQ = ez+£{B°{x)+£B^{£,x)) 
y p{{),x)/pQ=p + £[q°{x) + £q^{£,x)) 

where pis a constant, the functions and v^,B^,q^ are bounded in and where the 0-th order 

initial data verifies : 

I V±.v°(x)=0 (13.1) 

( 3/(z) such that B^{x) = f{z) — q°{x) (13-2) 
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then the solution of the full MHD system ([^ exists for a time T independe nt of e and this solution con¬ 
verges in to the solution of the reduced system given below in section 3.2.3 


Proof. The conditions on the structure of the system given in theorem 1 are satisfied while fhe 
conditions (13) express fhe facf fhaf fhe 0—order initial dafa is in fhe kernel of fhe large operafor. 
The assumptions of fheorem 1 are fhen fulfilled and fhe resulf follows. □ 


3.2.3 Slow limit of the system 

According to the general theory described in section 2.1 the solutions of (^ will be close fo fhe 
solufions of fhe limif sysfem of equations given by 


A^{Q)drW° + A,{W\ = 0 




(14) 


The zero-order solutions are functions = {v^, Vy,Vz,Bx,By,Bz, q) that are in the kernel of 
fhe large operafor. These functions musf fherefore verify : 


Vj.-u_L=0 (15.1) 

Vi(g + S,)=0 (15.2) 

using fhese resulfs, fhe explicif form of sysfem (14) can be written 

• V_l)u 2 ] -h d^q + (Sj^ • V±)q = 0 (16.1) 

^B^ + (u^ • Va_)B^ - {Bi_ ■ Vi)ui - d,vj_ = 0 (16.2) 


p {P)[-^vx + {vx ■ Vi)u^] VieV2 - d,Bj_ - {B^ ■ V^)B^ 


+ Vx(g' + 6i) 


= 0 


(16.3) 




B, + (ui • Vi)S, 


{Bi_ ■ Vj.)u2 -h V_L 


= 0 


(16.4) 


1 d 

+ {v±.^±)q] -P dzV^ -h 
or 


^±-v] 


= 0 


(16.5) 


where p and p{p) are consfant. Using fhe facf fhaf by (l^2) q + B^ = f{z) where f{z) is an 
arbitrary function, equations (1^4) and (1^5) can be combined to eliminate the corrective term 
Vl • v]_ resulting in : 


+ {B± • V_l)?;z + dzv^ = 0 

jp or 


We also note that in the perpendicular momentum equation, the term Vj^ (g^ 3- Bl) ensures that 
= 0, this term can therefore be combined wifh fhe V±B'^ 12 ferm wifh no change in fhe 
resulf. Infroducing fhe nofafions 


D-^ 

~Di' 


_d 


■ + {v±- V_l) • 


^//■ = {B±-y±)-+d,- 
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we get the final limit system : 


D - 

p—v_L-y //B± + V_lA =0 

(17.1) 

—Bj_ - V//v^ = 0 

(17.2) 

P^«z + V//g = 0 

(17.3) 

, 1 

(— - + = 0 

7P Dt " 

(17.4) 


where A stands here for a scalar "pressure" fhaf ensures fhaf fhe perpendicular divergence of 
fhe perpendicular velocify is zero. 

The equafions | [T7| shows fhaf fhe limif sysfem spiffs info fwo differenf sub-sysfems : 

• |[T7|1) and ([T7|2) as well as fhe consfrainf |[T5|l) describe fhe incompressible d 5 mamics of 
fhe perpendicular mofion of fhe plasma. This sef of equafion does nof depend on fhe 
pressure and foroidal velocify equafions and can be solved independenfly of fhe ofher 
fwo equafions. 

• On fhe ofher hand, fhe fwo scalar equafions 1^3) and (T7|4) describe fhe compressible 
parallel dynamics of fhe plasma. Acfually wifhouf fhe perpendicular convecfive ferms, 
fhese fwo equafions describe a compressible one dimensional flow in fhe parallel direcfion 
fo fhe magnefic field. Nofe fhaf fhese equafions are "slave" of fhe firsf fwo ones since 
bofh fhe perpendicular advecfion and fhe V// operafor depend only on fhe solution of 
equafions (T7| l) and |[T7|2). Thus |[T7|3) and |[T7|4) can be solved once fhe solufions of 
|[T7|1) and ||17[2) have been compufed. 

As in fhe original MHD sysfem, fhe sysfem (TT) is endowed wifh an involution : Using fhaf 
V.B = 0, we have in fhe limif e —0 fhaf fhe perpendicular divergence of fhe magnefic field is 
zero, 

Vj. • Sj. = 0 

if fhis properfy is frue for fhe inifial dafa, if is conserved by sysfem 

Proposition 1. Assume that the perpendicular divergence of the perpendicular magnetic field is zero at 
time f = 0 ; V = 0) = 0 then V±.B±{x, t) = Ofor t > 0. 

Proof. This follows direcfly by applying fhe perpendicular divergence operafor fo fhe perpen¬ 
dicular Faraday law ([T7|2). Nofe fhaf fo obfain fhis resulf, bofh fhe properfies V±.B± = 0 and 
= 0 are used. □ 

Alfhough, equafions (T7|l) and |[T7|2) have a similar sfrucfure, we nofe fhaf = 0 is not 

an involution for fhe sysfem : equafion |[T7|l) does nof conserve fhe perpendicular divergence 
of v±, fhe correcfive ferm Vj_A is fherefore needed fo insure fhaf Vj_ • = 0. 


We will now from fhe limif sysfem | [T7[ | obfain a reduced sysfem characferized by a smaller 
number of equafion fhan fhe number of fhe original sysfem. As explained in secfion 2.2 fhis 
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can be obtained by canceling out the corrective term. Since equations <[T7|l) and <[T7|2) form 
an autonomous system, we concentrate on these two equations. According to the general pro¬ 
cedure sketched in section 2.2 we look for a paramefrizafion of fhe funcfion space where fhe 
solufion belongs fo. In fhe presenf case, fhe space K = {(u_l, S_l); Vj_ • v± = • S_l = 0} can 

be paramefrized by 2 scalar funcfions cj), ijj such fhaf 


U_L = Gz X \7(j) (18.1) 

S_L = Oz X Vtjj (18.2) 

Lef us define for any scalar funcfion F G fhe operator A1 wifh values in x L^hy : 

M{F) = Gz X VF 


The following Green formula : 


G, X VF • Wdx = 


Fe„ X W ■ nds — / Fg, • V x Wdx 


Ion 


In 


shows fhaf fhe adjoinf operator of Ai is defined by : 

M*(W) =ez-V X W 


Using fhe general recipe given in secfion 2.2 we gef a reduced sysfem for fhe variables ((), ■0 by : 

D-^ 

pM*—M(0) - A4*V//M(0) = 0 (19.1) 


M*—MW-M*V//Mic/>)=0 (19.2) 

where fhe correcfive term V x A have been canceled ouf by fhe annhilafor operator A4 *. Affer 
some algebra, fhis sysfem admifs fhe following expression : 




( 20 . 1 ) 


drJ — V']_[dx(j)dy'ijj — dxipdycj)) — dJA = 0 (20.2) 

where U and J are defined asU = and J = note fhaf U — dyV^ — dxVy rep- 

resenf fhe z—componenf of fhe curl of fhe velocify vector, fherefore in reduced MHD liferafure, 
U is defined as fhe vorticity and |p0|l) is called fhe vorficify equafion by analogy wifh fhe fluid 
d 5 mamics case. 


From a physical poinf of view, fhe quanfify J = corresponds fo fhe toroidal currenf 

fraversing fhe plasma column and fherefore equafion <p0}2) defines fhe behavior of fhe foroidal 
currenf. In fhe framework of reduced MHD model, fhis equafion is nof used. Insfead, rewrifing 
@2) as 


- Vi + {dx4>dyi> - dxi>dy(j)) - dz(j)] = 0 (21) 

and nofing fhaf {dx4’dy'ijj — dxtpdycj)) corresponds fo fhe advecfion ferm ux • Vx0, one prefers fo 
use fhe equafion: 

—0 + ux-Vx0-^0 = O (22) 
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Strictly speaking < [2^ cannot be deduced directly from ( [2T) and integration factors should have 
appeared in However, it is possible to establish directly This is done in Annex 1. 


To complete the description of fhe reduced MHD models, we menfion fhaf in fhe presenf 
model, if is nof necessary fo solve fhe toroidal and pressure equafions p7[3) and p7|4) since fhe 
d 3 mamics is enfirely governed by p7|l) and p7|2). Neglecting fhese equafions, is also some¬ 
times justified as follows ( Il28l l : The acceleration ferm of fhe toroidal momenfum equation 
([^3) is : 

V //q = B± ■ V±q + -^q = B -Vq 

Then if can be shown (see 1281 ), fhaf if af fime f = 0, S • Vg = 0, fhen fhis quantify will sfay 
equal fo zero. Therefore, fhe foroidal accelerafion is null and if inifially Vz = 0, fhen fhe toroidal 
velocify will remain zero. Consequenfly, fhe velocify source V //Vz in fhe pressure equation 
remains zero and fhe pressure correcfion q behaves as a passive scalar. 

In fhe framework of MHD sfudies in fokamaks, fhe assumpfion B • Vg = 0 is very nafural since 
fhe flows under investigation are close fo an MHD equilibrium characterized by : 

Wp = JxB (23) 


fhaf implies fhaf B ■ Vp = 0. Acfually a lof of MHD sfudies aims fo examine fhe linear or 
non-linear sfabilify of such equlibrium and fherefore fhese works use precisely fhe relation 1231 
fo define fhe initial conditions. 

Summarizing, fhe d 5 mamics of fhe MHD model can be reduced fo a sysfem of 2 equafions 
for fhe scalar quantifies {(j), ip) 


P^U-ypJ = 0 (24.1) 

wifh 

u = -vl(p j = -v‘iip 

These equafions are conventionally written in a somewhaf differenf form emphazing fheir 
hamilfonian character 120] . Infroducing fhe brackef 


[f,9] = ez • Vj^/ X Vj.g 


we have fhaf for any / 

v± ■'^±f = [fj] while =[ipj] 

and fhe previous sysfem can be wriffen as 

- ^U + [f,ll]-[iP,J]-fj = 0 (25.1)' 

= 0 (25.2) 

where we have assumed p = 1 using an appropriafe scaling of fhe densify. 
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3.2.4 Fast modes of the system 


In this section, we briefly comment on the solutions of the full MHD system that are eliminated 

by the reduced model. In other term, we analyze the short time behavior of system |[T2). Con- 

* ' (2 

sidering the fast time scale r' = er or in an equivalent manner the fast reference time t' = —, 

i—1 

it is seen that the system |[12| reduces to the linear h 5 rperbolic system 


Ao{0)dr'W + Y,Cjd^,W = 0 

3 


(26) 


Let n = {rix, riyY be a 2D unit vector in the poloidal plane, the matrix ^o(O) + n-yCy) 

is diagonalizable and its eigenvalues are : 


Ao = 0 (with multiplicity 6), A+ = a / hP + 1 ^ nP + ^ 


(27) 


or in term of non-normalized variables : 


Ao = 0 (with multiplicity 6), A+ = 


HP' 




A_ =-i 


np- 


Bl 


(28) 


Comparing these expression to l|^, it is readily be seen that the non-zero eigensolutions corre¬ 
spond to fast magnetosonic waves traveling in the direction perpendicular to the toroidal mag¬ 
netic field BqBz. The situation here is quite similar to the one encountered with the compressible 
Euler equation where the fast limit corresponds to the acoustic equations describing the propa¬ 
gation of acoustic waves. Here, however, we also have an additional splitting in term of space 
directions. The fast limit of the system describes the propagation of magnetosonic waves in the 
poloidal plane while waves traveling in the toroidal direction are not present in this limit. The 
slow limit of the system that have been examined in section 3.2.3 thus excludes perpendicular 
magnetosonic waves in the same way as acoustic waves are filtered out from the compressible 
Euler equation when one consider the incompressible limit equation. 




Eigure 2: Comparison of the fast modes between the low Mach number limit and reduced MHD 
model; Left, Low Mach number limit : 3D isentropic propagation of acoustic waves; Right, 
reduced MHD models : 2D propagation of fast magnetosonic waves in the poloidal plane. 
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From a numerical point of view, this is one of the main advantage of reduced MHD since 
the use of the full MHD system |j^ implies strong CFL stability requirement linked to the prop¬ 
agation of magnetosonic waves. Note however that this splitting of the waves is not due to 
differences in the speed of propagation as for the Euler equation but rather to the different space 
scales in the toroidal and poloidal directions. In a toroidal system as a tokamak is, gradients in 
the toroidal directions are small with respect to gradient in the perpendicular directions and it 
is this fact that produce the wave separation rather than their speed of propagation since the 
velocities of Alfen and magnetosonic waves are roughly of the same order of magnitude. 


4 Concluding remarks 


This work has shown that the derivation of reduced MHD models for fusion plasma can be 
formulated in the general framework of the singular limit of h 5 rperbolic system of PDEs with 
large operator. This allows to use the results of this theory and to prove rigorously the validity 
of these approximations. In particular, it is proven, that the solutions of the full MHD system 


converge to the solutions of the reduced model displayed in section 3.2.3 


This work can be extended in several different directions. 


Eirst, the reduced MHD model considered in this paper is the simplest of a whole hierarchy 
of models of increasing complexity. The model used in the present work has at least two impor¬ 
tant weaknesses : 

a) It uses the straight tokamak model and therefore curvature terms are absent from the result¬ 
ing equations. More elaborated models 1^ 1121 retaining curvature effects and high order terms 
in e are available and can be possibly analyzed within the present framework. 

b) Another weakness of the model is that it uses as small parameter the ratio a/Ro that cannot 
be considered as small in a large number of today's machines. More elaborated models denoted 
in several references as "generalized reduced MHD models" iriTllT7ll30ll26l have been derived. 
These models do not make use of the small aspect ratio hypothesis and thus are in principle 
applicable with no restriction on the geometry. However, even from the point of view of formal 
asymptotics, these models are not always easy to understand and contains ad-hoc assumptions 
that are difficult to justify rigorously. It would be extremely valuable to study the possibility to 
formulate these "generalized reduced" MHD models along the lines exposed in this work. 

In the terminology of 1251 , the present work has examined the slow singular limit of the 
MHD equations. A second possible and interesting sequel of this work would be to examine 
the fast singular limit where no assumption is made on the boundedness of the initial time 
derivatives. On physical grounds, the assumption underlying the use of reduced MHD models 
is that fast transverse magnetosonic waves do not affect the d 5 mamics on the long time scale in 
the same way as in fluid d 5 mamics, the propagation of acoustic waves do not modify the aver¬ 
age incompressible background. Eor the Euler (or Navier-Stokes) equations this can be proven 
for certain cases e.g l25l in. Such a result however appears significantly more difficult to obtain 
for the MHD equations since their degree of non-linearity is higher than in the fluid d 5 mamics 
case. Note however, that the formal asymptotic expansion used in llTl can be considered as a 
first step in this direction. 
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A Annex 1 


In this section, we give a direct obtention of Equation (j2^. Since S is a divergence free vector 
field, fhere exisfs a vecfor pofenfial A such fhaf V y. A = B. From fhe expression ||^3) of fhe 
curl operator and fhe expression of fhe magnefic field, if is seen fhat ip is fhe toroidal component 
of fhis vecfor pofential. In term of vecfor pofential A, Faraday's law writes : 

d 

at 

where E is fhe elecfric field and f is fhe elecfric pofentiaQ 
Taking fhe scalar producf of fhis equafion by e^, one has 

^ip + e^- E- d;,(p = 0 

Now, using Ohm's law E + v x B = 0 and fhe identify 


—Bz ■ {v X B) = V ■ {e^ X B) 


one obfains (|22|| : 


d d 

—Ip + U_L • V_lV' - -7^(p = 0 

OT OZ 

Nofe fhaf since • W±ip = —B± ■ Vj_^, fhis equafion can also be wriffen 


From a physical poinf of view, fhis interprefafion shows fhaf fhe velocify defined by <[T8|l) is fhe 
so-called elecfric driff ve = E x B/\B^\ Indeed if can be shown (see ESil for insfance) fhaf fhe 
reduced MHD approximafion implies fhaf fhe fransverse elecfric field is elecfrosfatic : 


E± = "V ±(p 

from which one can deduce by faking fhe cross producf of Ohm's law by B fhe expression (1812) 
since in the small aspect ratio theory, the parallel and toroidal direction are identical up to terms 
of order e. 


^Note that the sign convention to define the electric field can be the opposite depending on the authors 
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